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SUMMARY 


A three-dimensional six-degree-of-freedom trajectory simulation applicable to 
preliminary stability and control studies of vehicles having large multi-engine boosters 
is formulated. The simulation employs rigid body equations of motion referenced to a 
special set of body axes known as Euler axes, an orthogonal right-handed triad having 
its origin at the instantaneous center of gravity of the vehicle. Fuel sloshing and elas- 
ticity effects are ignored. Simplifications are kept to a minimum in the derivation of 
the basic equations; consequently, the complexity of the simulation can be varied by 
assigning zero values to certain parameters. The aerodynamic characteristics of the 
vehicle are assumed to be a function of both Mach number and angle of attack. Vehicle 
mass and inertia characteristics are assumed to be time variant. Winds of arbitrary 
velocity and direction are incorporated into the simulation. 

Artificial stabilization is assumed to be provided by thrust vector control and 
application of the drift-minimum or load-minimum principle, although other control 
modes can be applied. Thrust forces are resolved in a manner which allows arbitrary 
location of the engines, except for the restrictions that the engine array must be sym- 
metrical with respect to the vehicle center line and gimbal points must lie in a plane 
perpendicular to the center line. 

The simulation will provide basic systems design information concerning the 
dynamic behavior of specific configurations with emphasis being placed primarily on 
the boost phase of flight. The basic information desired is the engine gimbal angles 
necessary for stabilization and the responsive motion of the vehicle to a specified 
disturbance. A numerical example is included which illustrates the application and 
flexibility of the simulation. Results are compared where possible with a typical two- 
dimensional simulation. 


SECTION I. INTRODUCTION 


Conventional body of revolution type vehicles are generally aerodynamically unstable 
since the center of pressure of the body alone is near the shoulder of the nose fairing- 
while the center of gravity is nearer the base. In the past, many vehicles have been*’ 
stabilized by the addition of fins to the body, or by a combination of fins and spinning 
the vehicle in flight (Ref. 1) . The provision of complete fin stabilization becomes less 
attractive with increasing size of the vehicle, for example those employed in manned 
space flight missions, because of large weight penalties and because the increasing fin 
span interferes with launching equipment. However, in some cases, relatively small 
fins may be used advantageously as a means of decreasing instability and reducing the 
load imposed on the main control system. Spin stabilization is rarely even considered 
for large vehicles, especially those using liquid propellants because of the resulting- 
fuel sloshing and centrifugal force problems. 

Artificial stabilization by means of a control system, rather than the provision 
of a stable airframe, appears to be the best approach for large vehicles (Ref. 2) . Such 
a control system requires control torques which may be produced by air vanes, jet vanes, 
engine gimbaling, and several other means or combinations of these. Control by engine 
gimbaling, sometimes referred to as thrust vector control, has been successfully em- 
ployed in several large vehicles and appears to be attractive for those of the future. 
Consequently, only gimbal control is considered in this report. The limiting of maximum 
required control deflections to reasonably small values is important in view of the in- 
creasing complexity of gimbal bearings, adequate closure of the vehicle's base for 
protection from base heating, and structural requirements for engine mounts and air- 
frame (Ref. 2). 

Since the success of a mission often depends on proper attitude stabilization of the 
vehicle, control requirements and restrictions need to be considered during the prelim- 
inary design phase when the basic vehicle configuration is established. Such studies 
should be as simplified as the desired degree of accuracy allows. Furthermore the 

methods used should be quite flexible to cope with the variety of configurations usually 
encountered. J 

A unique method of artificial stabilization of aerodynamically unstable vehicles 
known as the "drift minimum principle" (Ref. 3) has been widely used in preliminary 
control investigations and is of considerable current interest (Refs. 4 and 5). Appli- 
cation of this scheme renders the vehicle path almost insensitive to wind forces, which 
are the primary source of disturbance during the early stages of flight. This principle 
utilizes linear control equations in which the so-called control gains are predetermined 
in a specified manner for each vehicle. A variation of the principle, in which temporary 
path instability is allowed in order to reduce the loads on the vehicle, is referred to as 
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the "load-minimum principle. " Both the load-minimum and drift-minimum control 
modes are discussed more recently in Reference 6, with two other possible control 
schemes. A more generalized discussion of the problems of attitude stabilization for 
large guided missiles can be found in Reference 2. 

Preliminary vehicle design applications of the above principles have been almost 
entirely restricted to two-dimensional, planar analyses. Consequently, these studies 
are limited to independent investigations of pitch, yaw, and roll motions* with roll 
motion frequently neglected entirely. Analyses of this type are partially justified for 
symmetrical configurations. For non-symmetrical vehicles, which could be the case 
for configurations carrying winged payloads for manned space flight, a three-dimensional 
analysis of the dynamic motion is necessary to study the aerodynamic and coupling effects 
properly (Ref. 7). In a three-dimensional simulation, pitch, yaw, and roll motions are 
examined simultaneously. For multi-engine boosters, three dimensions are often neces- 
sary to resolve the thrust forces and moments properly. 

In the two-dimensional analyses, flight conditions must be assumed constant in all 
planes of motion except the plane being considered. Often the vehicle aerodynamic 
characteristics are assumed to be constant (Ref. 2) although the Mach number and angle 
of attack may be changing. Furthermore, wind disturbances in the two-dimensional 
analyses are necessarily restricted to constant direction and it is assumed that a direct 
side wind results in the maximum yaw disturbance. In Reference 8, however, it was 
found that in most cases, the maximum 3 ^aw dispersion does not coincide with a direct 
side wind. Winds of arbitrary velocity and direction are easily incorporated in a three- 
dimensional simulation (Ref. 9). 

Restricted three-dimensional equations, with linear thrust vector control, are 
formulated in Reference 10, but are not suitable for studies involving multi-engine 
vehicles. Previous studies involving multi-engine vehicles have been limited to a 
particular engine arrangement (Ref. 4). The literature reviewed indicates that there 
is a need for a three-dimensional trajectory simulation, applicable to stability and 
control studies, incorporating variable wind, variable aerodynamics, and arbitrary 
engine arrangement. 

This study formulates a general three-dimensional six-degree-of-freedom trajec- 
tory simulation applicable to stability and control studies of vehicles having large multi- 
engine boosters. The simulation is formulated in a manner allowing arbitrary location 
of the engines with respect to the center line of the vehicle, except for the restrictions 
that the engine array must be symmetrical, and gimbal points must lie in a plane per- 
pendicular to the center line. The vehicle is assumed to be a rigid body relying on 
artificial stabilization by thrust vector control. Elasticity and fuel sloshing effects 
are ignored. 
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The simulation will provide basic systems design information concerning the 
dynamic behavior of specific configurations with emphasis primarily on the boost phase 
of flight. The basic information desired is the engine gimbal angles necessary for 
stabilization and the responsive motion of the vehicle to a specified disturbance. A 
linear control system is assumed and two control modes, namely, the drift-minimum 
and load-minimum principles (Ref. 6), are considered. 

A ficticious, but realistic, vehicle is postulated as an illustrative example and 
typical' results are presented. The results are compared, where possible, with typical 
two-dimensional results. Indications are that the simulation derived in this study 
satisfies the intended purpose and is adequate for preliminary stability and control 
studies of large multi-engine vehicles. 


SECTION II. COORDINATE SYSTEMS 


A. COORDINATE DESCRIPTION 

The choice of coordinate systems and the methods of coordinate resolution 
are extremely important in the formulation and solution of any system of equations. 
Several coordinate systems used in solving the three-dimensional flight equations have 
been investigated in Reference 11. The conclusion is reached that the best axis system 
from the standpoint of simplicity is one employing wind axes for translational equations 
and body axes for rotational equations. Such a procedure is rarely employed, however, 
because confusion is inevitable when different axis systems are used for the two sets of 
equations. 

Wind axes have long been used in analytical stability and control studies (Ref. 12) 
to take advantage of the resulting simplifications in the equations of motion. However, 
since wind axes are not fixed with respect to the body, it is necessary either to assume 
that the time rate of change of the moments and products of inertia is negligible, 
formulate an auxiliary expression for their time dependency; consequently, a special 
set of body axes known as Euler axes is used in this study. 

The Euler axis system (Ref. 13) is associated with an orthogonal right-handed 
triad having its origin at the instantaneous center of gravity of the vehicle as shown in 
Figure 1. 

The X-direction is positive along the longitudinal axis pointing forward. The Y 
and Z directions are the pitch and yaw axes, respectively, with the Y-axis positive to 
the right looking in the flight direction. With respect to the Euler axis system, the 
vehicle has linear and angular velocities and accelerations but no displacements. This 
system is especially adaptable to studies of rigid body flight dynamics since the velocities 
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FIGURE 1. ILLUSTRATION OF VEHICLE BODY AXES 
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measured with respect to these axes are the same as the velocities that are measured 
by instruments mounted in the vehicle. Another advantage of this system, as opposed 
to wind axes, is that the moments and products of inertia are independent of time, 
except for changes in the vehicle's mass. 

Another set of axes used in establishing the equations of motion is referred to 
as earth-fixed axes. A flat, non-rotating earth is assumed for this study so that the 
earth-fixed axis system is actually an inertial system which is needed when applying 
Newton's Laws. For a flat earth, these axes are at all times parallel to the axes of the 
vehicle's stabilized platform. 

In the earth-fixed system (Xg, Yg, Zg) , the Zg axis is the vertical axis being 
positive in the opposite direction of gravity. The Xg axis is defined by the desired 
flight azimuth with the Yg axis completing the right-handed orthogonal system. 

A third system, often needed in connection with aerodynamic parameters, is 
the stability or experimental axis system. This system is commonly used in wind 
tunnel facilities and the like for referencing aerodynamic data. However, depending 
upon model instrumentation, such data can also be referenced to body axes. For this 
study, it is assumed that all aerodynamic parameters are referenced to the body or 
Euler axis system. 

B. COORDINATE RESOLUTION 

Before going into a derivation of the equations of motion, it is desirable to 
develop the transformations relating the body and earth-fixed coordinate systems. As 
yet, no parameters have been defined for relating the two. For this purpose, the Euler 
angles 9 , , and [ are introduced and will be referred to as the pitch, yaw, and roll 

angles, respectively. Euler devised a method of rotating an axis system successively 
about its own axes so that at the end of the third rotation the system is parallel to a 
second specified axis system. Care must be taken in using Euler angles since they do 
not obey the commutative law thereby requiring that a specified sequence of rotation be 
maintained throughout a particular analysis. 

A standardized Euler sequence has been used in formulating the three-dimensional 
flight equations for aircraft and unguided vehicles (Refs. 7, 9, and 11). This particular 
sequence has the order yaw, pitch, and roll with the pitch angle referenced to the hori- 
zontal plane. By a twist of fate, all Euler sequences have a singularity when the second 
rotation is exactly +90° . The singularity is a zero in the denominator of one of the 
equations relating the angular velocity components of the vehicle to the time rate of 
change of the Euler angles (Ref. 11). Consequently, the conventional sequence used 
for studies involving aircraft and unguided vehicles has a singularity at 9= +90° , and 
is obviously unsuitable for vertically launched vehicles. As a result of this situation, 
a different Euler sequence is used in the present study. The rotational order for this 
study is 0, , and : (pitch, yaw, roll) with 0 being referenced to the vertical Zg axis. 
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It is pointed out in the development of the transformation equations that this system has a 
a singularity at ^ = +90° , a condition that is not likely to occur. 

An arbitrary orientation of body and earth-fixed axes is shown in Figure 2 with 
the Euler angles indicated. Also shown are the axes about which the time rates of 
change of the Euler angles are measured. It should be noted that these axes are not 
mutually orthogonal, thus explaining why the Euler sequence is not commutable. In 
their original positions, X coincides with Zg and Z coincides with Xg while Y and Yg 
are 180° out of phase. Any arbitrary orientation of the body and earth-fixed axes is 
obtained as follows: 

1. Rotate about the Yg-axis through the angle 9 to define the X ls Y 1} Z t system 
(Fig. 3). 

2. Rotate about through the angle f to form the X 2 , Y 2 , Z 2 system. 

3. Rotate about Z 2 through the angle $ to establish X 3> Y 3 , Z 3 . The X 3 , Y 3 , and 
Z 3 axes correspond to the body axes Z, -Y, and X, respectively. 

The angles are positive when, in performing the rotations described, the right- 
hand rule is applied. 


From Figure 2, keeping in mind that initially Y and Yg are 180* out of phase, it 


can be seen that: 


Xi = 

x g 

COS 

6 - Z„ sin 9 

O 

Yi = 

Y g 



Zi = 

z g 

cos 

Q 4 - Xg sin 0 

X 2 = 

Xi 



y 2 = 

Yi 

cos 

^ + Zi sin i f 

Z 2 = 

Zi 

cos 

- Yi sin \|r 

x 3 = 

z = 

: x 2 

cos $ + Y 2 sin 0 

y 3 = 

-Y 

= Y 2 cos 3> - sin $ 

z 3 = 

X = 

z 2 



> (la, b, c) 


> (2a, b, c) 


> (3a, b, c) 
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9 cos 4 j COS (j5 T 



Step 3 

X 3 , Y 3 , Z 3 Correspond to Z,-Y, X Respectively 


FIGURE 3. DETAILS OF EULER SEQUENCE 



Combining Equations 1, 2, and 3 gives 


X = Xg sin 0 cos t " Yg sin \|r + Zg cos 0 cos \]/ 

Y = Xg (cos 0 sin $ - sin 0 sin \|r cos 0) - Yg cos cos $ 

- Zg (cos 0 cos <I> sin + sin 0 sin $) )>-(4a, b, c) 

Z = Xg (cos 0 cos <t> + sin 0 sin \| r sin $) + Yg cos \|r sin $ 

+ Zg (cos 0 sin \]/ sin <p - sin 0 cos $) 

Equations 4a, 4b, and 4c represent the geometric resolution of body and earth-fixed 
axes. In matrix notation (Ref. 15), the transformation is 


(5) 




\ 


\ 

X ] 


dn ^12 di3 | 


fxg 

Y 

I 

= 

d 2i d 22 d 23 


Yg 

z i 


1 ^31 d3 2 d33 j 


1 Zg / 


where 

dn = sin 0 cos \|r 
di 2 - - sin ^ 
di 3 = cos 0 cos \|r 

d^i = cos 0 sin 0 - sin 0 sin cos 0 
d 2 2 “ - cos cos <2> 


d 2 3 = - (cos 0 cos <!> sin \|r + sin 0 sin $) 
d3i = cos 0 cos $ + sin 0 sin i|f sin 0 
d32 = cos \|r sin 0 

d 33 = COS 0 sin f sin <£> - sin 0 cos 0 


These quantities are actually the direction cosines of angles between the earth-fixed and 
body axes expressed in terms of Euler angles. The methods for deriving the direction 
cosine relations by means of spherical trigonometry are outlined in Reference 7, but the 
approach used above presents fewer difficulties. Reference 7 also gives the orthogonality 
relations 
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d±i 2 + d± 2 ^ + d X 3 2 = 1 

^21^ + d22 2 + d 2 y ; = 1 

^ 31 ^ + d 32 “ + d.33'- = 1 


y 


(6a, b, c) 


These relations provide the means for a partial check of solutions of the equations of 
motion since the values determined for the Euler angles must satisfy the equations. 

The matrix operation given in Equation 5 can be inverted to transform quantities 
from body to earth -fixed axes. Thus, 


Xg^ 


/dn 

^21 

X g 

= 

di2 

& 22 . 

Z8 < 


^di 3 



d3i \ 

d32 

d 3 3 


(7) 


Tne geometric transformation can only be used when transforming translational 
quantities such as distance, velocity, and acceleration. In the development of the 
equations of motion, the time rate of change of the Euler angles will be of interest. 
Another look at Figures 1 and 2 shows that the actual axes about which the rotations take 
place are Z 2 , Yg, and X t . The rotational resolution of this triad and the body axes is 
a necessity. The angular velocities about the X, Y, and Z axes are denoted by p, q, and 
r, respectively, and represent three rotational degrees-of-freedom of a rigid body. The 
relative orientation of the angular velocity vectors in the two systems is illustrated in 
Figure 2. From the figure, it is evident that 


p = <t>-0sin\lr 

q=\|rsinO - 0 cos cos $ 
r = cos <t> + 0 cos sin 


y 


(8a, b, c) 


J 


These equations represent the rotational resolution of the triad X, Yg, Zj to the 
body fixed axes. In matrix notation, 


(9) 


1 p) 


h 

- sin \|r 

o \ 


q 

= 

0 

- cos cos <& 

sin $ 


r j 


1° 

cos \|f sin 0 

COS <$> j 



Inverting this matrix gives 


ii 



I 0 = 


(r sin <t> - q cos cl)) 


cos \]f 

\| f = r cos 3> + q sin 0 


V (10a, b, c) 


I ^ = p + (r sin 0 - q cos $) tan \|/ 
or i 


/ . \ 

\ 


\ 


* 1 


1 - cos <t> tan \|r sin $ tan \[r | 


p j 


e 

= 

q cos 0 sin <f> 


q 

(id 



cos \|r cos a|t 




Wi 


^0 sin 0 cos 0 j 


i r i 



Equation 10a has a singularity point at = +90°. This substantiates the earlier 
discussion pointing out that every Euler sequence has a singularity when the second 
rotation is exactly +90° . No further transformations are necessary for this study. 


SECTION m. EQUATIONS OF MOTION 


A. BASIC EQUATIONS 

The approach used in establishing the equations of motion is first to formulate 
the general case under certain restrictions, and consider simplifications later. For 
this study, the vehicle is considered a rigid body in space. Gravity is assumed to be 
constant in both magnitude and direction. The earth’s rotation and orbital motion are 
neglected, making the earth-fixed axes an inertial system. 

The assumption that the earth is fixed in an inertial reference frame ignores the 
existence of the Coriolis and centripetal accelerations due to angular velocity of the 
earth. Actually these two components of acceleration usually comprise only a small 
portion of the total acceleration experienced by the vehicle. Since the inertial frame 
is attached to the earth, additional errors result as flight progresses over the surface 
of the earth. Since the assumed gravitational field is in error by the inclination of the 
gravity force vector, neglecting the curvature of the earth's surface will introduce an 
altitude error. If the flight time is not excessive, and flight occurs over a small part 
of the earth's surface, all these errors may be neglected. These conditions are not 
met by either ballistic or complete orbital trajectories. For preliminary dynamic 
stability and control investigations of the boost phase of large vehicles, the area of 
interest is normally limited to atmospheric flight and the assumptions above appear 
justified. 
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The quantities u, v, and w are actually the speeds of the center of gravity. The 
shift of the center of gravity during flight, with respect to some fixed vehicle reference, 
will be accounted for in the equations of motion; however, the time rate of change of 
the center of gravity (the velocity of the center of gravity with respect to the vehicle) 
is small and can be neglected. 

Newton's laws state that the sum of the external forces and moments acting on 
a rigid body can be equated to the time rate of change of linear and angular momentum, 
respectively. Accordingly, the force and moment equations are: 




+ k ; F, 


( 12 ) 


dt 



j 


M + 


k ) N 


(13) 


where Ly and Hy are the linear and angular momentum, respectively, F x , Fy, F z , 
and L, M, N are the force and moment components along and about the X, Y, and Z 
axes, respectively. 

The subscript T in the equation above refers to the total system and indicates 
the inclusion of gases in the nozzle generated by burning of propellant. To take into 
account the dynamic effect of the jet on the vehicle, the rate of change of linear momen- 
tum must include the rate of linear momentum leaving the nozzle exit in the jet. Likewise, 
the rate of change of angular momentum must include the rate of change due to gases 
leaving the nozzle. 

The rigid body equations of motion for flight mechanics applications under all the 
assumptions made earlier, but with the additional restriction of constant mass, are 
derived in References 10 and 13. Vehicle mass is variable during the boost phase of 
any vehicle. The effect of mass loss has been considered by many investigators in 
analyses concerning such problems as the performance of high altitude sounding rockets 
and escape from the earth by a rocket (Refs. 14 and 15). In these analyses, however, 
dynamic stability was not included in the area of investigation. 

It has been determined, however, that the variation in mass causes so-called jet 
damping forces and moments. A restricted two-dimensional analysis has been made in 
Reference 15 to determine the effect of mass variation on the dynamic stability of rocket- 
propelled vehicles. The vehicle involved in that study was a small, low-thrust, single- 
engine rocket. The results show that a disturbance of the variable mass vehicle damps 
out more rapidly than it does with the mass assumed constant, indicating that the variable 
mass vehicle is more stable. The effect of jet-damping is dependent on the ratio of the 
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time rate of change of mass to the gross mass of the vehicle. Since large multi-engine 
vehicles are generally aerodynamically unstable, any stabilizing effect should certainly 
be considered as a means of alleviating the control problem; therefore, mass was 
assumed variable in this study. 


The basic three-dimensional six-degree-of-freedom equations are derived in 
Reference 16 for a single-engine missile. It is shown that Equations 12 and 13 can be 
expanded into the form 

dL dL _ _ 

~ = dt + * (V " V e + Q X r e } < 14) 

and 


dt 


dH _ _ _ _ 

— + mr x (V + ft x r ) 
dt e e e 


(15) 


where V e is the velocity of the exhaust gases relative to the vehicle, V is the trans- 
lational velocity of the vehicle, Q, is the rotational velocity of the vehicle, and F e is 
the vector distance from the. instantaneous center of gravity to the center of the nozzle 
exit. The mass flow rate, m, is defined as 

. dm 

m = ' dt (16) 


Furthermore, L and H are the linear and angular momentum neglecting the exhaust 
gases. 

Following the methods of Reference 17, remembering that L and H are referred 
to the body axes system which is rotating with an angular velocity ~ with respect to 
the earth-fixed or inertial system, the scalar components of Equations 14 and 15 for the 
case of a multi-engine vehicle can be shown to be 
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where the subscript n refers to a particular engine. 

These equations represent the motion along and about the X, Y, and Z axes. The 
equations are complete, within the assumptions made to date, except for the external 
forces and moments which appear on the left side of each equation. 


One immediate simplifying assumption can be made concerning r e ^ where, as 


defined above, 


n 


ir + i r +kr 
ex ey ez 

n n 
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For single engine boosters where the nozzle center line coincides with the longi- 
tudinal axis, r e y and r ez are identically zero. For multi-engine vehicles with the engines 
mounted symmetrically about the longitudinal axis. 


n n 



Actually Equation 23 is ti’ue only when the engines are not gimbaled. Neglecting the 
small bias resulting from gimbaling or engine-out conditions, and assuming that the 
engines always lie to the rear of the center of gravity, then 

r = -i r 
e ex 

n n 


If it is further assumed that all nozzle exits lie in a plane approximately parallel to the 
Y-Z plane, then r ex is the same for all engines and the subscript n can be dropped. 
Under these assumptions, the basic equations reduce to 
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( 29 ) 
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Further simplifying assumptions could be made at this point. For example, it 
could be assumed that the X - Z plane is a plane of symmetry so that Iyz = Ixy = °* 

This condition is normally a basic design criteria; however, in some cases, it may be 
desirable to determine the effect of certain symmetries. In studies where analytical 
or closed-form solutions are a desired result, it is necessary to linearize the equations 
(Ref. 18) . Such situations are restricted to a solution for a particular type of problem 
and are not readily adaptable to a variety of unexpected problems which typically arise 
in the preliminary design phase. Where large highspeed electronic computers are 
available, the best approach is to hold simplification to a minimum; then if simplications 
are desired, certain parameters are simply set equal to zero. In cases where the 
simplified equations cannot handle a particular problem, the zeros are easily removed 
and no major modifications are necessary. Consequently, no further simplifying 
assumptions are made concerning Equations 24 through 29. The purpose of this study 
is not to establish any particular stability criteria (Ref. 19) , but to establish a flexible 
simulation applicable to the variety of problems encountered in preliminary design. 

The complete equations of motion assumed in this study are established by formu- 
lating the various forces and moments represented by the left-hand side of Equations 24 
through 29. For this purpose, it is assumed that 


i ; F x + 3/ f y + V. F Z ' F a + F T + F G + F M 


i)L + OM+k)N-M M + M A + M T 


(30) 

(31) 


where 


F M , m m ~ miscellaneous forces and moments such as those resulting from thrust 
misalignment or retro-rockets 

F q = gravitational force 

Fa» M a = aerodynamic forces and moments 

Fy, M-p= forces and moments caused by thrusting action. 
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It should be noted that there are no gravitational moments, since the origin of the 
body axes is at the instantaneous center of gravity of the vehicle. Each of the forces 
and moments will now be developed separately according to their source. 

B. GRAVITY FORCES 

The gravitational vector is referenced to the earth-fixed system and is 
directed along the negative Zg axis. The transformation given in Equation 5 can be 
used to resolve the vehicle weight into body-axes components as follows: 


F = - mg cos 9 cos 

G x 

F = mg (cos 0 sin ^ cos 0 + sin 0 sin 0) 
G 

Y 

F = mg (sin 9 cos $ - cos 9 sin sin $) 


(32) 


where g is the magnitude of the acceleration of gravity. 

The introduction of 9 , , and 0 into the equations of motion adds three unknowns 

to the six basic unknowns defined earlier in Equations 24 through 29. The situation is 
such that these are six equations with nine unknowns. Consequently, auxiliary equations 
are needed to determine 9 , f , and <t> in terms of the basic unknowns. In Equation 10 
it was shown that the time rate of change of the Euler angles can be expressed in terms 
of the rotational velocity components p, q, and r. The Euler angles can thus be cal- 
culated from the relations 

9 = / 9 dt + 9 

o 

f = / f dt + f 

o 

d> = / $ dt + $ 

O 

where Q q , t Q , and $ 0 are initial values and 9 , ¥, and f> are given in Equations 10a, 

b, c. 


y (33a, b, c) 
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C. AERODYNAMIC FORCES AND MOMENTS 


The aerodynamic forces and moments can be formulated in terms of the 
normal and axial force coefficients and the so-called stability derivatives. The aero- 
dynamic forces and stability derivatives are usually referenced to either body or 
experimental axis system, depending on model instrumentation or methods of calculation 
(Refs. 9 and 13). In this study, it is assumed that all aerodynamic data are referenced 
to the body axis system since the basic equations of motion have been established in 
that system. Following the approach of References 7, 12, and 13, it is further assumed 
that all stability derivatives for large missiles, with or without non-symmetrical payloads, 
are either zero or negligible except 


Axial Force 



Side Force 



Normal Force 



Rolling Moment 



C 


z 


r 


Pitching Moment 



, C 


m 


q 


Yawing Moment - c , C 
n- n 
a y r 

where a P and ay are the pitch angle of attack and yaw angle of attack ( commonly referred 
to as the angle of sideslip) as illustrated in Figure 4. These angles will be discussed 
later. 


The stability derivatives and the experimental and analytical methods for deter- 
mining them are discussed thoroughly in the literature,' particularly in many reports 
published by the National Advisory Committee for Aeronautics (NACA) and later by the 
National Aeronautics and Space Administration ( NASA) . The methods of calculating 
the stability derivatives will not be discussed here; however, in any stability investigation 
it is desirable to have specific analytical or empirical relations for calculating these when 
experimental data are not available. 

In more detailed stability and control studies, or when an attempt is made to 
establish a general stability criteria (Ref. 19) , it may be necessary to account for more 
of the numerous stability derivatives. For this study, those listed above are sufficient. 
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Standard NASA notation is employed ( see List of Symbols) in defining the stability 
derivatives. The so-called static derivatives are simply partial derivatives with respect 
to a particular variable. For example 


where is the normal force coefficient. 

The dynamic stability derivatives, in which the variable in the denominator of the 
derivative is associated with motion, and nondimensionalized by multiplying the partial 
derivative by the parameter 2 Vr/D where Vr is the resultant flow velocity. Thus 


m 

q 



dc 

m 



where C m is the pitching moment coefficient. 

The convention is adopted that drag, side, and normal forces are positive along 
the negative X, Y, and Z axes. Thus the aerodynamic equations are 
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where CP and CG denote the distance from a reference plane passing through the engine 
gimbal points to the center of pressure and center of gravity, respectively. Furthermore 
S and D are the reference area and diameter, respectively, and the dynamic pressure, 
q', is calculated from 



where is the free-stream density and is the resultant flow velocity. 

The aerodynamic forces and moments are seen to be dependent upon the resultant 
air-flow velocity. In the absence of wind the velocity of the vehicle relative to the 
earth-fixed system, V", and its velocity relative to the flow, , are equal in jnagnitude. 
Since the basic equations of motion are derived in terms of the components of V, it is 
necessary to calculate Vr to account for the possibility of wind. 


Wind data are usually referenced to the earth-fixed coordinate system. The wind 
velocity can be broken into components X^, Y W g, and Z W g along the Xg, Yg, and Zg 
axes, respectively. These components can be resolved into the body axes system, using 
Equation 5, and combined algebraically with u, v, and w to yield the respective components 
of Vr as follows: 
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c) 


In the presentation of wind data it is customary to list vertical gusts, Z W g, and 
the resultant wind parallel to the Xg - Yg plane, necessitating the calculation of X^g 
and ^ W g. A scheme similar to that presented in Reference 12 is convenient for wind 
resolution in the Xg - Yg plane. The resultant wind parallel to this plane is denoted 
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by Vw, whichis referenced to true north by the direction angle ( Fig. 5) . The Xg 
axis, corresponding to the desired launch azimuth, is referenced to north by the azimuth 
angle, A. 

From Figure 5 it is obvious that 


y v ■ A 


w 


where 7^ is introduced as a means of referencing to the Xg axis. Now it follows 
that 
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wg 
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(36a, b) 


The introduction of winds during launch creates large angles of attack (+90° for 
vertically launched vehicles) which are beyond the range of linear aerodynamics. Non- 
linear aerodynamic characteristics with respect to angle of attack are assumed to 
simulate the launch motion under wind conditions. Furthermore, since Mach number 
varies considerably during the flight of a vehicle, it is also necessary to consider the 
variation of the aerodynamic characteristics with Mach number. Mach number is 
determined by the usual relation 

V R 

M = — (37) 

V S 

where Vg is the speed of sound. 

Referring again to Figure 4, the pitch and yaw angles of attack are seen to be 
defined by 
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(38a, b) 
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The resultant angle of attack is determined from 



where u^, v^, and are components of the resultant flow vector. 

The time derivatives of the pitch and yaw angles of attack are needed in the 
calculation of the aerodynamic forces and moments. From Equations 38a, b: 
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cos 2 OL 
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(40a) 
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D. THRUST FORCES AND MOMENTS 

Before deriving the equations for thrust forces and moments, it is necessary 
to adopt a sign convention and define several angles associated with the engine. For a 
multi-engine booster, a number system must also be adopted to facilitate references to 
a particular engine. For this purpose consider the typical engine arrangement for an 
eight-engine booster shown in Figure 6. 

In such an arrangement, the outer engines are usually the control engines while 
the inner ones are held fixed. No restriction has been placed on the relative locations 
of the control and fixed engines. For reference purposes, the control engines were 

assigned odd numbers (i, 3, 5, 7, ) and the fixed engines were assigned even 

numbers (2, 4, 6, 8, ). The only other restriction is that the complete arrange- 

ment will be symmetrical with respect to the center line of the vehicle. 

The angles associated with the engines, in both the initial and gimbaled positions, 
are defined in the list of symbols. Using engine one in Figure 6 as an example, the 
angles are illustrated in Figure 7. Note that the initial position illustrated is typical for 
all engines, whereas the gimbaled position is applicable only to the control engines. 

The sign convention adopted is such that positive Py-p an( * ^Pt pro^ 1106 force 
components in the positive Y and Z directions, respectively. 
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The thrust force and moment equations can be derived by using engine one of 
Figure 7 as an example, then amplifying to obtain the general equations for a multi- 
engine vehicle. 


From Figure 7 note that 


Pv ~ Py + 

T x 01 Y 1 




(41) 

(42) 


where Py^ an d are the yaw and pitch gimbal angles measured from the original 

canted axis of the engine. Note that an initial cant angle is not a restriction, but is 
considered as the general case since rocket engines are quite often canted. The main 
purpose of this is to direct the thrust of the engine toward the center of gravity at 
liftoff, thus reducing moments that could occur in the event of engine failure in a multi- 
engine vehicle. 

Further relations which are either obvious from Figure 7 or easily derived are 
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Gimbaled Position 


FIGURE 7. ILLUSTRATION OF ANGLES ASSOCIATED WITH INITIAL AND 
GIMBALED POSITIONS OF ENGINE NUMBER ONE 



The angle P n is introduced as a means of quickly locating the thrust direction with 

respect to the center line of the vehicle. Looking along the X axis in Figure 7b, |3p^ 

is the angle between the X-Z plane and the plane containing the thrust vector and the X 
axis. From this figure it is seen that 
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The component thrust forces are easily developed for the engine in its initial 
position, being simply 
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where Tj is the thrust of the example engine. 


The development of the forces for the engine in a gimbaled position is more 
difficult. If the forces are formulated in their exact form, it is difficult to study the 
resulting motion independently in pitch and yaw because of the dependency of P ^ on 
p and • For example, the exact form of the Fp^ i n a gimbaled position is 
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The dependency of Fp^ on ^ Yi anc * ^Pf obvious. As will be seen later in 
the development of the control equations, it is desirable to have the engine pitch force 


dependent only on the pitch gimbal angle £ , and known initial values. The same 

applies for having the yaw force component dependent only on known initial values and 


the yaw gimbal angle, £ y. • The desired relationship for Fp can be obtained by 

l Zj 

assuming that its dependency on P y^ is negligible for reasonably small angles. For 
this first approximation, Equation 54 becomes 



The desired result hqs been obtained since Fp^ is now dependent only upon the initial 

values and (3 . The dependency is of a trigonometric form, however, making it 

difficult to see readily the change in Fp because of a change in . 

Zf 1 

Suppose then, as a second approximation, it can be assumed that Fp can be 

Z 1 

represented by its value in the engine's undeflected position, given by Equation 55, plus 
a linear gradient through P y 1 = Pp A = 0. In other words, it is assumed that 
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and then shown that the last term can be neglected within certain limits so that finally 
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The gradient of Fy with respect to Pp^ is obtained by partial differentiation of 
Equation 54 with the result 
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For reasonably small values of the angles, this gradient can be assumed linear. At the 
initial point of Py^ = Pp = 0 the value is 
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Making use of Equations 44 and 47, this is further reduced to 
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For small initial angles, the cross coupling term represented by Equation 58 can be 
ignored in comparison with the related term given by Equation 57. As a result, Equation 
56 can be applied so that 
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This expression of F .p is also dependent only upon initial values and the pitch gimbal 




angle, Pp^. I n this case, as opposed to Equation 55, the dependency of F^ on 
is a simple linear relationship since the initial values are known and the coefficient of 
P is a constant. The relative accuracies of both approximations must be compared 
with the exact value given by Equation 54. To make this comparison, typical values for 

the initial cant angle must be assumed. Consider the case ^Yqi = ^Pol = 3 °’ Knowing 
these two values, the remaining initial values can be computed. Figure 8 shows the 
comparison of the approximate and exact values over a wide spread of values for /S Yi . 
As indicated in the figure, f3 values of five and ten degrees are assumed which 
correspond to the lower and upper set of curves, respectively. The difference between 
the first and second approximations, representing Equations 55 and 59, is seen to be 
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FIGURE 8. COMP ARIS ION OF EXACT AND APPROXIMATE METHODS FOR CALCULATING F T 




















less than one per cent for Pp^ = 10° , with the first approximation being closer to the 
exact value. Because of this small difference, the second approximation is used to 
take advantage of the straightforward expression. 

Referring again to Figure 8, it can be seen that for values of Py^ U P t° ten 

degrees, the error involved in using the second approximation amounts to only three 

per cent for Pp^ = 10° , and less than three per cent for Pp^ = 5°. For values of 

Py^ exceeding ten degrees, the error increases at a moderate rate. A gimbal 

angle of ten degrees is conservative for most vehicles and use of the second approx - 
imation is justified. 

The remaining forces and moments caused by thrust are developed under the same 
assumptions. Following the procedure applied previously to the pitch force, the yaw 
force is found to be 
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The exact form for the axial component of thrust is 



For this component, the pitch and yaw gimbal angles share equal importance so that the 
one cannot be neglected in favor of the other. To avoid this double dependency, both 
P and Pp^ can be neglected for the angular magnitudes normally encountered. The 

axial component thus becomes 
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which, by use of Equation 45, reduces to 


34 


= T cos R 
1 P 01 


( 61 ) 


F 


1 


Equation 61 is identical with Equation 51 which expressed the value of in the engine 

undeflected position. Ignoring the effect of Py^ anc ^ u P on F<j« * s tantamount 

to assuming that the change is negligible in the axial component of thrust because of 
gimbaling. 


The development of the thrust force component for the example engine has been 
completed. It should be emphasized that, although a particular engine is used in the 
development, the results are applicable to any engine no matter what its location relative 
to the vehicel center line may be, so long as the sign convention specified is used. Thus 
for the general case Equations 59, 60, and 61 can be written as 
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where the subscript n refers to a particular engine (n = 1, 2, 3, 4, ) . Earlier it 

was specified that all fixed engines would be identified by even numbers. Therefore for 
n, an even number, the gimbal angles Py and Pp n are zero and terms involving these 
variables have no meaning. The summation of components for all engines of a multi- 
engine vehicle can be written 
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where n again refers to any particular engine, control or fixed, and n c (n c = 1, 3, 5, ) 

refers to control engines only. Note that the first summation term in both Equations 63 
and 64 will be zero if the engine arrangement is symmetrical with respect to the center 
line and all engines perform as expected. In case of engine failure or thrust misalignment, 
however, biases are created which result in an unbalance of forces and moments. Thrust 
misalignments are not considered in the development of the thrust equations but a term is 
included in the general equations of motion to account for such misalignments. 


Since all the initial values are known in Equations 62 through 64, the equations can 
be expressed in a simpler notation by grouping certain constants. For this purpose, let: 
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Then Equations 62 through 64 can be written in the simplified form 
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The components of the moment caused by thrust are developed in a manner similar 
to that used above for the force components. In the development of the pitching moment, 
for example, it is assumed that 


d M x 
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ii Pi P] 


^ die re engine number one is again used for illustration. The gradient term is again 
as.sumed linear and evaluated at the initial condition = Furthermore, 

cro53S-coupling terms are found to be small for the gimbal angle range of interest and 
cons equently are neglected. Thus pitching moments caused by thrust are assumed to 
be dependent only on pitch gimbal angles and yaw moments dependent only on yaw gimbal 
anglers. It is also assumed that, except for an engine malfunction, the summation of 
thrust moments are zero when none of the engines are in a gimbaled position. Under 
these assumptions, the pitch and yaw components are found to be 
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where R n is the perpendicular distance from the vehicle center line to the gimbal point 
of a particular engine, and 
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4n 

The production of rolling moments with thrust vector control is also necessary 
for asymmetric configurations to counteract induced aerodynamic rolling moments. 

In a later development of the control equations, it is assumed that when only pitch forces 
and moments are produced by thrust gimbaling, all of the pitch gimbal angles of all 
control engines are equal. The same applies to yaw forces and moments produced by 
yaw gimbal angles. If this restriction is maintained in the overall control of the vehicle, 
it is impossible to produce rolling moments, except in the case of a malfunction, because 
of the earlier restriction concerning the symmetry of the engines. Therefore to produce 
coordinated rolling moments, certain pairs of engines must gimbal with equal magnitudes 
but in opposite directions. This will be taken into account in the control equations. For 
practical multi-engine configurations, roll rates should not be severe so that large 
rolling moment corrections are not needed. Thus it can be assumed that roll control 
can be accomplished by either pitch or yaw gimbaling independently. 

Following previous procedures, the general expression for thrust rolling moment 
was found to be 
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where and C 2 are so-called "go-or-no-go" constants. For pitch gimbal roll control, 
Cj = 0 and C 2 = i and vice-versa for yaw gimbal roll control. Also 

K = cos 8 tan sin (3 

Y on P p 


39 



“l 

2TI, 

3n 

4n. 


cos 8 sin 6 
on K p 

c on 


cos 2 p 


1 - sin 2 p 


on 


on 


cos 8 tan 8 cos 8 
on Y p 

on on 



cos 8 cos 8 
on p 

c on 


cos" p^ 


on 


1 - sin 2 p sin 2 6 
on p 

c on 


The thrust in all previous equations is assumed to be a function of time to simulate 
thrust build-up and decay with altitude. The altitude correction is given by the usual 
relation 
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where T on is the thrust of a particular engine at some reference pressure P Q , A en is 
the nozzle exit area, and P w is the trajectory free stream static pressure. It is 
assumed that the thrust of all engines is measured at the same reference pressure, 
and that each functions ideally or not at all. 

E. MISCELLANEOUS FORCES AND MOMENTS 

In the formulation of the general equations of motion, it is virtually impossible 
to account for all possible sources of external forces and moments. Forces and moments 
caused by thrust misalignments, fin deflections, or retro-rockets, as a few examples, 
are not considered in detail. Consequently, as a means of making the equations as general 
as possible, a miscellaneous term is added to each equation. Thus, 
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The miscellaneous terms can be either constant or time variant. Actually the term 
F X]VI is needed in the axial force equation to balance the gravitational force at time zero. 

At this point, as will be seen later, the axial force equation reduces to 

F x M = mg (72) 

F. GENERAL EQUATIONS OF MOTION 


The general equations of motion are obtained by summing the X, Y, and Z 
components of the various forces and moments, and substituting on the left-hand side 
of Equations 24 through 30. The general equations are 
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= pl x + qr (I z - y + (pr - q) 1^ - (pq + r) 1^ + (r 2 - q 2 ) I yz 
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Table 1 is included as a convenient reference, listing all inputs to the equations 
according to the independent variables of which they are a function. 
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TABLE 1 


PROGRAM INPUTS 


Inputs which are a function of time only are 

F X m ’ F Y m > F Z m ’ L M’ M M’ N M’ m ’ m ’ I X 5 I Y’ I Z’ I X’ V I Z 

I XY’ I YZ’ I XZ’ I XY’ *XZ’ T on’ A en’ CG ’ r ex’ 3 0Y 

a 0P’ a lY’ S 1P’ b 0Y’ b 0P’ b 0Y’ h 0P’ ®c’ Q c’ 


Inputs which are a function of Mach number only are 


Inputs which are a function of altitude only are 


V , V , Z 
s w wg 


Inputs which are a function of Mach number and pitch angle of attack, or its derivative, 
are 

c , c , c Y , c , c 

m m Z m. 

^p 

Inputs which are a function of Mach number and yaw angle of attack, or its derivative, 


c v c v c v c \' \ 


Inputs which are a function of total angle of attack are 
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TABLE i (CONCLUDED) 
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Initial values of the basic variables are 
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SECTION IV. CONTROL EQUATIONS 


The six general equations of motion representing the six degrees of freedom of a 
rigid body given by Equations 73 through 78 contain 11 unknown variables, namely the 
linear velocities (u, v, w), the rotational velocities (p, q, r) , the orientation angles 
( 9, a , 0 ) and the gimbal angles ( 3 y , Pp ) • The auxiliary relations, given by 

n c 

Equations 33 a, b, and c, determine the orientation angles. Two additional equations are 
needed to match the number of equations with the number of unknowns. These additional 
relations are the control equations which can be employed for the determination of the 
gimbal angles. 

The control equations have a tremendous influence on the eventual motion of the 
vehicle. The physical behavior of a vehicle without control varies considerably during 
its flight because of variation in the vehicle's mass and its aerodynamic characteristics. 
The control system has to be adapted to these variations which are often not accurately 
known. Control requirements for a given configuration are dependent on the mode of 
operation of the control system. This will be discussed in more detail later. 

One of the basic assumptions of this study is that the vehicle is considered a rigid 
body. However, deformation of the vehicle under the influence of aerodynamic and 
maneuver loads, together with propellant sloshing, creates interactions between the 
associated modes of motion which require more complex analysis. Interaction between 
motions around the three body axes is usually small since the control system keeps 
angular velocities and deviations small. For preliminary design purposes the control 
equations can be assumed to contain only linear functions of certain control variables, 
with coupling in the control system ignored. Consequently, it was assumed that yawing 
motions are controlled by yaw gimbaling and pitching motion of pitch gimbaling. As 
mentioned earlier, it is assumed that roll control can be accomplished by either pitch 
or yaw gimbaling. In this study, it is further assumed that the vehicle is equipped with 
an attitude reference system and an angle -of-attack meter so that the control variables 
are the attitude angles and angle of attack. Under these assumptions, the general control 
equations are 
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Path control, in which displacement variables and their derivatives are included 
among the control variables, usually has little or no effect on the gimbal angle require- 
ments and can be neglected in initial studies ( Ref. 4) . 

The control equations are idealized by neglecting derivatives of the gimbal angles 
on the left side which produce phase lags but are not important for basic studies (Ref. 2). 
Furthermore, all derivatives of the control variables are neglected in this study except 
the first derivatives of f and 9 . The control equations can now be written in this form 

n 

c 
n 

c 

where the coefficients of the control variables are the so-called control gains and the 
variables with the subscript c indicate the desired value of that variable when the desired 
value is not zero. If the desired value of a variable is assumed to be zero at all times, 
no such term is necessary. As an illustration, consider the yaw control equation. If 
a solution of the general equations of motion indicates the vehicle is yawing, the yaw 
angle f is amplified by the gain factor aQy and a yaw gimbal angle is created to correct 
the yawing motion. In the pitch plane, however, the trajectory is normally tilted so that 
the desired value of the pitch angle, Q c , is not zero at all times. Thus if a solution of 
the equations of motion should indicate that the actual pitch angle e is not equal to e c , 
then a pitch gimbal angle is initiated and the trajectory will be corrected. The desired 
value of the roll angle, $ c , is normally zero except in cases where launch equipment 
restricts the alignment of the vehicle so that a programmed roll is necessary. All the 
control gains and the desired trajectory parameters are time variables. 

Using Equations 79 a and b for control, the simulation established in this study 
can be used for two specific types of problems. In the first case, completely arbitrary 
control gains can be assumed and the resulting dynamics of the vehicle determined. In 
the second case, which is the one encountered more often, the control gains are established 
according to some preset relation determined by the desired motion of the vehicle during 
flight. The establishment of a control program by which an otherwise unstable vehicle 
is stabilized is commonly referred to as artificial stabilization. As indicated earlier, 
the control modes capturing the most interest at present are the drift-minimum principle 
( DMP) and the load-minimum principle ( LMP) . Both principles are discussed in detail 
in References 3 and 6. 

The DMP is simply a control program which minimizes the dispersion of the vehicle, 
in both the pitch and yaw planes, in response to a wind shear. This is accomplished by 
forcing the vehicle to assume an attitude which leads to the approximate cancellation of 


‘ “OY * + \ * + V °Y + h 0Y 
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forces perpendicular to the undisturbed flight path. The wind induced oscillations and 
accelerations approach zero and a component of thrust against the wind cancels the lift 
and side forces. Although a detailed derivation of the DMP is given in References 3 
and 6, a short derivation is included here because of differences in sign convention and 
basic definitions. In the development of the DMP, it is assumed that the vehicle is not 
rolling and that control in the pitch and yaw planes can be analyzed independently. Con- 
sidering then the motion in the pitch plane with a wind shear imposed on the vehicle, it 
is assumed that the wind is horizontal so that the resultant wind in the pitch plane, for 
small values of ? r , is X w g ( Fig. 9) . 


Referring to Figure 9, let 

9-9 = A0 
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where s c and Xg c are the desired pitch angle and the desired rate of change of range, 
respectively. Also from the figure, for small angles 
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where a 
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is the wind angle, having the approximate value 
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Consequently Equation 81 can be approximated by 
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The incremental motion in the Xg direction can be approximated by 
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FIGURE 9. ATTITUDES OF AERODYNAMICALLY UNSTABLE VEHICLE IN 

RESPONSE TO WIND WITH ATTITUDE AND ANGLE OF ATTACK 
CONTROL 
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In the derivation of the DMP, it is assumed that rolling motion is small and can be 
neglected; hence Py an d Pp are the same for all engines. 

The incremental pitching motion can be approximated by 
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so that Equations 83 and 84 can be written in the form 
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From Equation 79 b 
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where rolling motion is again neglected. 

Substituting Equations 82 and 87 into Equation 86 gives 
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Substitution of Equations 82, 87, and 88 into Equation 85 gives 
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For DMP it is assumed that 
Xg = AS = A9 ~ 0 
so that Equation 89 becomes 
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Equation 90 represents the drift minimum condition and is identical with that derived 
in References 3, 4, and 6 except for certain differences in signs. These sign differences 
are solely caused by the sign conventions used here and do not detract from the basic 
principle involved. Equation 90 can also be derived as shown in Reference 6 by establishing 
the characteristic equation of Equations 86, 87, and 88, and examining the roots in detail 
as a check on the validity of the assumptions made. Such detail is not necessary here, 
since the basic principles are well known. 
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The DMP condition, represented by Equation 90, is not sufficient for determining 
the magnitude of the gains a Qp an d bQp but only indicates the linear relationship between 
them. Consequently another relation is needed to determine their values. By virtue of 
the DMP condition the translational path dispersion is practically eliminated (Refs. 3 and 
6) so that the rotary motion of the vehicle about its center of gravity is predominant. 
Equation 89 is a linear differential equation describing the rotary motion in response 
to wind disturbances. From this equation, the undamped natural frequency of the rotary 
motion is found to be 
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The frequency has a definite influence on the control requirements. Care must be 
taken to insure that the frequency is not so low that it corresponds to the bending 
frequencies of the vehicle to give resonance effects. The high values of frequency are 
bounded by the limitations of the control system hardware. Frequencies currently being 
used range from about 0. 2 to 0. 9 cycles - per - second. As is pointed out in Reference 2 
the optimum control frequency for a particular vehicle depends on the characteristics of 
the wind disturbances encountered. Once the frequency is established, or assumed, 
Equations 90 and 92 can be combined to determine the values of aQp and bQp. Adding 
the two equations gives 
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Once a G p is determined Equation 92 can be used to calculate bop. 

The value of a^p for the DMP is determined by the damping of the rotary motion 
about the center of gravity. The ratio of damping to critical damping is found from 
Equation 88 to be 
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that if the frequency and damping are known a iP is readily determined. 
*vf f is 0. 75. 


so 

value of £ 
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The lateral control gains for the DMP are established by identical procedures with 
those employed above. Thus analagous to Equations 91, 92, 93, and 94, it was foun 
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When the DMP is applied, attitude deviations are small and the resultant angle of 
attack does not approach zero when wind disturbances are encountered. As a result, 
normal loads are imposed on the vehicle. In some cases these loads may be beyond 
the design limit of the vehicle and the DMP cannot be employed. The angle of attack 
can be reduced by reducing a oy and a 0P since this allows greater path instability. The 
extreme is reached when a OY and a OP are equal to zero (Ref. 3) . In this situation, the 
vehicle turns its nose in the resultant flow direction like a weathercock stable vehicle 
without special attitude control so that attitude deviations of the vehicle may be large 
Since the angle of attack, and thus the aerodynamic loads, approaches zero int is ex reme 
case, this is commonly referred to as the load-minimum principle ( LMP) . This approach 
may be indicated for large, highly unstable vehicle configurations where tremendous 
control torques and structural loads are involved ( Ref. 2) . If proper guidance equipment 
is available, the deviation of the flight path in the presence of wind is of little concern 
since it can be compensated for by guidance induced maneuvers during the latter part o 
propelled flight when the vehicle is relatively insensitive to aerodynamic effects. 
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OY are set equal to zero in Equations 92 and 96, it is found that 
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Note that the load-minimum principle is not a zero load 
the least amount of load that can be associated with the DMP. 

however, there is a specific frequency at which the LMP does 
condition. 


condition but represents 
For a given vehicle, 
reduce to a zero-load 


, • e s r ong influence of C lp , C 2P , C lY , and C 2 y on both the DMP and LMP is 

obvious. The parameters C 1P and C 2P are the specific aerodynamic restoring torque 

nd the specific control torque (pitch plane) , respectively, while C 1Y and C 2Y are 

identical parameters for the yaw plane. The fact that all four parameters are determined 

olav heti 6 ° haracter * stics of the vehicle and its trajectory is illustrative of the inter- 
play between design and control. 


Since it is assumed in the development of the DMP that roll can be neglected the 

T r™'r‘ h °* and h °P Also no guarantee haaleen rnide 

that the DMP or LMP gains satisfy the stability conditions of the servomechanisms 

employed m the control system (Refs. 20 and 21). This requires a more detailed 
anaiysis. In view of this uncertainty it is assumed that arbitrary values for the roll 
g ins can be chosen in any particular investigation, with the actual values being deter- 
mined by the limitations of the control system hardware. 

... J enera l equations of motion are formulated so that roll can be controlled by 

ither pitch or yaw gimbaling. In some cases it may not be necessary to use all engines 
for roll control. For this reason the control equations are now written in the form 
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C > + h Qp ($ - 0 ) (102) 


where the constants K R ^ are roll proportionality constants taking on values of zero, or 

phis or minus one. If the value is zero then the engine to which it is assigned does not 
participate ss the roll control. Once it is established which engines will be used for roll 
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. t hp analyzed to determine which will be assigned 
control, the engine or one value simply causes certain engine 

a plus or a minus one value. T p . . . nnnos ite directions so that roll 

pairs to deflect with equal incremental magnitude but in opposite direct. 

restoring moments are initiated. 


SECTION V. SOLUTION OF THE EQUATIONS 


Once the vehicle data is established and the toittal 
general Equations 73 through 78 can ' Je ^^earth-fixed acceleration components 

s -7 vr body 

axis components should include both linear and rotational accelerations. 
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The earth-fixed velocity and trajectory geometry components are obtained by 
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SECTION VI. NUMERICAL EXAMPLE 
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parameters, for the drift-minimum and load-minimum control modes. Typical result 

vehicled and c ° mpared > where Possible, with two-dimensional results for the same 
hide characteristics. Past experience has shown (Refs. 4 and 51 that th* • 

STa bal 3 tf 6 requirements occur when a vehicle is disturbed by a high velocity wh^T™ 
r , m ® ^S lon °f maximum dynamic pressure. Consequently the numerical 

vehicfe e t rTjectory! S ‘ ri ° ted * ““ ” aXta “ m PreSSUre re * ion ° f the burned 

A vehicle having eight engines arranged and numbered as shown in Pirnre fi t„ 

n<i7 U i me ’-f W1 thC f ° Ur outboard engines having gimbal capability. The thrust and 
nozzle exit area of all engines are the same, and the gimta! points arainTpl tT 

1 500 000^™ H VeMclC , Center line ' 1116 total sea-level thrust of the vehicle is 

thrust build- up" is assnmTf mg t ° 1 ' lft ° ff We ‘ ght ° f 1 ' 000 ' 000 P^nds. An instantaneous 
duration Of 120 1! rf * or su »P llc «y so that T on is constant with time. A thrust 
duration of 120 seconds is also assumed. The initial cant ancle 0 arp 4° an H r 

the inboard and outboard engines, respectively. From Figure 6, the ’following angles^ 
defining the initial thrust direction of each engine are obvious. 


P = 45° 

p 

P n = 

225° 

01 

P 

05 


0 

0 

Ch 

II 

Q 

QCl 


270° 

02 

06 


0 = 315° 

P os 

6 

‘ P 

07 

135° 

B = 0 

11 

c 

ca 

180° 


The outboard engines lie on a circle about the vehicle center line such that 
Rl = R 3 = R 5 = r 7 = 8, 5 feet 
Similarly, 

R2 = R4 = Rg = Rg = 3* 5 feet 

(Eo 7 R «‘! C °” tr01 WU1 be provided by y»w deflections of the control engines so that 
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From Figure 6 note that a yaw deflection of engines one and three in the negative 
direction would tend to correct a positive roll, while a positive yaw deflection would 
be needed for engines five and seven. Thus, assuming that all of the control engines 
share in roll control, the roll proportionality constants are (Eqs. 102 and 103): 


3 

-*» - 1 

7 

The mass, center of gravity, inertia, and thrust moment arm characteristics of 
the assumed vehicle are shown in Figures 10 and 11. It is assumed that all the mis- 
cellaneous forces, such as Fy^jj are zero. The assumed aerodynamic characteristics 
of the ficticious vehicle are shown in Figures 12 through 15. Since the current example 
is limited to the region of maximum dynamic pressure, the angles of attack will not be 
large and most of the curves shown do not consider the variation of the aerodynamic 
parameters with angle of attack. Figure 13 illustrates the more general case, although 
the curves shown for large angles of attack are not needed for this example. A double 
iteration procedure is employed in the simulation when both the Mach number and angle- 
of-attack variation of the aerodynamic parameters are considered (Ref. 9) . 

The desired (no wind) trajectory parameters in the region of maximum dynamic 
pressure are given in Table 2. The desired roll angle, : c , is zero for the entire flight. 

The desired flight path angle, Q c , and its first time derivative are the only parameters 
shown in Table 2 that are required inputs to the control equations. 

The desired launch azimuth is due east from Cape Canaveral, Florida so that 

A. = 90° 

Typical wind profiles and wind shear data for Cape Canaveral are given in Reference 
22. A typical wind profile is shown in Figure 16. It is assumed that the vehicle formulated 
in the present examti y . must be designed to fly through the wind shear represented by this 
profile. The assumed wind profile is designed so that the maximum wind velocity corre- 
sponds to the same altitude at which the vehicle encounters maximum dynamic pressure 
(Table 2) . The wind shear is assumed to have an 8. 5 second build-up; consequently, 
the vehicle enters the shear at t = 51. 5 seconds, which is assumed to be the starting 
time for calculations made in this example. The initial conditions of the desired tra- 
jectory are thus 


K R. 
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Flight Time, t ( seconds) 


FIGURE 10. MASS, CEN TER OF GRAVITY AND 
THRUST MOMENT ARM DATA 
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Mass, m (slugs) 
















FIGURE 11. MASS MOMENTS AND PRODUCT OF 
INERTIA CHARACTERISTICS 
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Roll Moment Of Inertia, lx (lb-ft-sec 2 ) 















Mach Number, M 

FIGURE 12. CENTER OF PRESSURE AND AXIAL FORCE COEFFICIENT 
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(dimensionless) 
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C Y 7 

y *y 

(1/rad) 
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( 1 / rad) 
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3 



Mach Number, M 

FIGURE 13. RATE OF CHANGE OF NORMAL FORCE 

COEFFICIENTS WITH ANGLE OF ATTACK 
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m *p 

(1/rad) 


Pa Y 
(1/rad) 
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( 1 / rad) 
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(1/rad) 



FIGURE 14. PITCHING MOMENT AND YAWING MOMENT 
STABILITY DERIVATIVES 
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FIGURE 15. ROLLING MOMENT STABILITY DERIVATIVES 
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Altitude, Z„ (feet) 



Wind Velocity, Vyy (ft/sec) 

FIGURE 16. WIND VELOCITY PROFILE 
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= 51.5 sec 


t 

o 

u = 1190.8 ft/sec 
o 

v = 0 
o 

w = 14.8 ft/sec 


X = 6920. 6 ft 
go 

Y =0 
go 

Z = 26191.0 ft 
go 


p 0 = 0 

=- 0.0091 rad/sec 

r = 0 
o 

f = 0 

o 

0 = 23.76° 

o 

0=0 

o 


The drift “minimum attitude and angle-of-attack control gains can be calculated 
from Equations 92, 93, 96, and 97, and are shown in Figures 17 and 18 for several 
undamped natural frequencies of the control system. The attitude control gains are 
zero for the load-minimum principle. The load minimum angle-of-attack control gains 
can be calculated from Equations 99 and 100, and are shown in Figure 19. The attitude 
rate control gains are the same for the LMP and the DMP but are a function of the 
damping ratio. In this example, it is assumed that the pitch and yaw damping ratios 
are equal. The attitude rate control gains calculated from Equations 94 and 98 are 
shown in Figures 20a, b, and c for several assumed damping ratios. 

Initial calculations are made for both drift-minimum and load-minimum control 
under the following assumptions: 


f 


nY 


nP 


0.5 cps 


£ 


Y 


0.75 


0 = 180° 

w 
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Flight Time, t (seconds) 


a 0Y 




OP 


FIGURE 17. DRIFT MINIMUM YAW AND PITCH ATTITUDE 
CONTROL GAINS 
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b 0Y 

b 0P 



Flight Time, t (seconds) 

FIGURE 18. DRIFT -MINIMUM YAW AND PITCH ANGLE OF 
ATTACK CONTROL GAINS 
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Angle Of Attack Control Gains, bgy And bQp (dimensionless) 


FIGURE 


b 0Y 



Flight Time, t (seconds) 


LOAD-MINIMUM YAW AND PITCH ANGLE OF 
ATTACK CONTROL GAINS 








Attitude Rate Control Gains, Attitude Rate Control Gains 
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FIGURE Z0. YAW AND PITCH ATTITUDE RATE CONTROL GAINS 
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Attitude Rate Control Gains 
a, Tr And a,„ (seconds) 



Flight Time, t (seconds) 


(c) Damping Ratios £,y And £,p Equal 1.0 


FIGURE 20. (CONCLUDED) 







The assumptions concerning the roll control gains h OY and h Q p are completely 
arbitrary. A direct side wind ( S w = 180°) is selected to facilitate the comparison of 
the current calculations with the results of a typical two-dimensional simulation. The 
two-dimensional analysis accounts only for normal and axial aerodynamic forces which 
are assumed constant. 

The drift-minimum control requirements are examined first. Figure 21 shows the 
time history of the gimbal angle and lateral acceleration in response to the assumed 
wind disturbance. In this and succeeding figures, time is assumed to be zero at the 
instant the vehicle encounters the wind disturbance. In other words, zero time corre- 
sponds to the actual flight time of t = 51. 5 seconds with maximum dynamic pressure 

occurring 8. 5 seconds later. The finite gimbal angle shown in Figure 21 for the present 
simulation at zero time is the value required to fly the desired trajectory without wind 
conditions. Calculations are stopped when the variables approach steady-state conditions 
Time histories of the yaw angle and angle of attack, for drift-minimum control are shown 
in Figure 22. 

Application of the load-minimum principle yields similar results as shown in 
Figures 23 and 24. Load-minimum control reduces the gimbal angle and angle of attack 
as compared with drift-minimum control. The yaw angle and lateral acceleration are 
larger and continue to increase rather than approaching a steady-state condition, thus 
indicating path instability. 

The results of the present analysis are in general agreement with two-dimensional 
results for both control principles. This indicates that the two-dimensional analysis gives 
satisfactory results when the disturbed motion of the vehicle takes place primarily in 
one plane as it does in the case of a direct side wind. 

The time histories of the gimbal angle calculated from the present simulation are 
not associated with any particular one of the four control engines of the assumed vehicles; 
rather the maximum value at each instant of time is used. However, since the gimbal 
angle relative to the initial canted axis is shown, the value is the same for all four engines 
except for small differences necessary for roll control. 

The roll angle of the assumed vehicle is negligible with both load-minimum and 
drift-minimum control for the case of a direct side wind. The dispersions of the flight 
path angle and longitudinal acceleration are shown in Figure 25. For the case of a side- 
wind the two-dimensional analysis neglects motion in the pitch plane; hence no comparison 
can be made for the pitch plane variables. 

All of the results shown thus far are for a frequency of 0. 5 cycles per second. The 
effect of frequency on the maximum values of lateral acceleration, yaw angle, gimbal 
angle, and angle of attack are shown in Figure 26 for drift-minimum control. Similar 
results are shown in Figure 27 for load minimum control. Values for yaw angle and 
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Present Analysis 
Two-Dimensional Results 



Time After Wind Disturbance (seconds) 


FIGURE 21. TIME HISTORY OF GIMBAL ANGLE AND LATERAL 
ACCELERATION IN RESPONSE TO WIND DISTURB- 
ANCE (DRIFT-MINIMUM PRINCIPLE) 
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Present Analysis 



FIGURE 23. TIME HISTORY OF LATERAL ACCELERATION AND GIMBAL 
ANGLE IN RESPONSE TO WIND DISTURBANCE 
(LOAD-MINIMUM PRINCIPLE) 
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Angle Of Attack, 

a (degrees) Yaw Angle, (degrees) 


Present Analysis 
Two-Dimensional Results 



Time After Wind Disturbance (seconds) 


FIGURE 24. TIME HISTORY OF YAW ANGLE AND ANGLE OF ATTACK 
IN RESPONSE TO WIND DISTURBANCE 

(LOAD-MINIMUM PRINCIPLE) 
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Dispersion Of Dispersion Of 

Longitudinal Flight Path 

Acceleration, Anlge, A0 

AX„ (ft/sec 2 ) (degrees) 


— Drift Minimum Principle 

Load Minimum Principle 



Time After Wind Disturbance (degrees) 


FIGURE 25. TIME HISTORY OF DISPERSION OF LONGITUDINAL 
ACCELERATION AND FLIGHT PATH ANGLE 


77 












Xl 

a 

< 


■Present Analysis 
Two-Dimensional Results 
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Natural Frequency, f n (cycles/seconds) 
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FIGURE 26. EFFECT OF NATURAL FREQUENCY ON MAXIMUM VALUES 
OF LATERAL ACCELERATION, YAW ANGLE, ANGLE OF 
ATTACK AND GIMBAL ANGLE (DRIFT-MINIMUM PRINCIPLE) 
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lateial acceleration are not shown for load— minimum control since they increase con- 
tinuously and consequently have no maximum value ( Figs. 23 and 24) . The two- 
dimensional results here are generally conservative as compared with the results of the 
present analysis. A notable exception is the yaw angle with drift-minimum control, 
the two-dimensional analysis predicting 20 per cent lower values ( Fig. 23) . 

The effect of wind direction on the lateral -acceleration is shown in Figure 28 for 
both control principles. Similar variations for gimbal angle and angle of attack are 
shown in Figure 29. The case of a failure of engine number one ( T t = 0) is also 
considered in Figure 29. All of the curves shown in Figures 28 and 29 have mirror 
images for the other 180-degree range of wind direction except for the curves con- 
sidering engine failure. 

It is apparent from Figures 28 and 29 that the maximum points on the curves are 
not coincident with a direct side wind ( 6> w = 180°) as is sometimes assumed. For the 
present example, however, the maximum values are not significantly different from 
those corresponding to a side wind. 

A significant roll angle results in the case of engine failure, for some wind 
directions. The range of the roll angle for 9 W between 90 degrees and 270 degrees for 
the assumed vehicle is as follows: 

- 1.20° < $ < - 2.76° 

It is beyond the scope of this example to consider all of the possible variations 
and combinations of the parameters involved. A complete stability and control inves- 
tigation of any vehicle requires a tremendous volume of charts and calculations. While 
the results shown are brief, they illustrate the application and flexibility of the simulation 
developed in this report. The simulation should be restricted primarily to first stage 
flight because of the assumptions of a non-rotating, flat earth and a constant acceleration 
of gravity. 
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FIGURE 29. EFFECT OF WIND DIRECTION ON MAXIMUM VALUES OF 
GIMBAL ANGLE AND ANGLE OF ATTACK 
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